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Standard physics is still capable to interpret ~18 TeV photons from GRB 221009A 
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ABSTRACT 


It is reported that the Large High Altitude Air Shower Observatory (LHAASO) observed thousands 
of very-high-energy photons up to ~18 TeV from GRB 221009A. We study the survival rate of these 
photons by considering the fact that they are absorbed by the extragalactic background light. By 
performing a set of 10° Monte-Carlo simulations, we explore the parameter space allowed by current 
observations and find that the probability of predicting that LHAASO observes at least one photons 
of 18 TeV from GRB 221009A within 2000 seconds is 4-5%. Hence, it is still possible for the standard 
physics to interpret LHAASO’s observation in the energy range of several TeV. Our method can be 
straightforwardly generalized to study more data sets of LHAASO and other experiments in the future. 


1. INTRODUCTION 


More than 5000 photons above 0.5 TeV emitted from GRB 221009A at redshift z ~ 0.1505 were observed by the 
Large High Altitude Air Shower Observatory (LHAASO) ! within 2000 seconds after the first detection by Swift, 
Fermi-GBM, Fermi-LAT, and so on °. The highest energetic photons were reported to reach ~ 18 TeV, representing 
the first observation of photons above 10 TeV from gamma-ray bursts (GRBs). Such observations intrigued studies 
on photon mixing with axion-like particles Galanti et al. (2022), Lorentz symmetry violation Li & Ma (2022), and 
both Baktash et al. (2022). In our work, we will investigate the survival probability of multi-TeV photons from GRB 
221009A by considering the fact that they are significantly absorbed by the extragalactic background light intervening 
between the GRB and the Earth. We will further show whether the standard physics is still capable to interpret 
current observations. 

Very-high-energy photons can dissipate their energies via annihilation with photons in cosmic microwave background 
(CMB) and extragalactic background light (EBL), producing electron-positron pairs. The threshold of this channel 
to happen is Ey, = m2/E,, where me and F, are the mass of electrons and the averaged energy of background light, 
respectively. Therefore, for CMB photons, this threshold is hundreds of TeV, implying that we can safely disregard 
the effect of CMB photons. However, the energy of EBL photons can be higher by several orders of magnitude than 
the CMB photons, changing the threshold to be lower by the same magnitude. For example, the threshold is ~ 2.6 
TeV if we consider Ey ~ 0.leV. Therefore, we should take into account the suppression effect of EBL photons on the 
detected flux of ~ 18 TeV photons by LHAASO. 


2. FLUX OF TEV PHOTONS AND EBL ATTENUATION 
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Figure 1. Posterior probability distribution functions of A; estimated in the case of priors a; = —1.87 + 0.04 (blue) and 


—2.12 + 0.11 (red), respectively. 


The EBL-suppressed flux F, depends on the intrinsic flux F; and the optical depth e~7 due to absorption by EBL. 
Therefore, we have 
P(E) = (Be (1) 
where F is the observed energy of photons and z is the redshift of GRB 221009A. We use the tabulated data of the 
EBL and optical depth measured by Ref. Dominguez et al. (2011). The intrinsic flux of photons is approximated to 
be a power-law MAGIC Collaboration et al. (2019) 


F,(E) = A; (x) (2) 


where A; and qai, respectively, stand for the intrinsic flux and spectral index at the pivot energy scale 0.5 TeV. 
Considering the reports of Fermi-LAT, we have two measured values of a;, namely, —1.87 + 0.04 ? and —2.12+0.11 
4 which were obtained at 0.1-1 GeV, but in different temporal intervals. We leave A; to be determined by LHAASO 
as follows. 

By considering the performance of LHAASO Cao et al. (2019), we predict the number of events with energy above 
0.5TeV to be 


Nso.stev = r | F(E)Sea(£, O)dE , (3) 
0.5TeV 


where Seg is the effective area of LHAASO-WCDA, as provided in Fig. 6 of Chapter 3 in Ref. Cao et al. (2019), 
0 ~ 28° is the zenith angle of GRB 221009A, and T = 2000 seconds is the duration of LHAASO observation. 

To explore the parameter space, we perform Bayesian analysis by considering the fact that the number of photons 
above 0.5TeV is more than 5000, as reported by LHAASO °. We assume that the event number follows Poisson 
distribution with probability distribution function (PDF), i.e. 


p(k) = \¥e~*/k! (4) 


with À = 5000 being the expectation and k = N>o.5Tey. Our results and conclusion are remained when choosing other 
value of A within [5000, 6000). We implement the Bayesian inference by using the affine-invariant Markov chain Monte 
Carlo (MCMC) ensemble sampler in emcee Foreman-Mackey et al. (2013). We assume that A; has a uniform prior in 
the range [1 x 10-4, 7 x 1074] in units of TeV-'m~?s~!. We take the aforementioned Fermi-LAT measurements of the 


3 https: //gen.gsfc.nasa.gov/gcn3/32658.gcn3 
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Figure 2. Probability of predicting that LHAASO-WCDA observes at least one photon from GRB 221009A within 2000 
seconds. The coloring of curves is consistent with that of Fig. 1. 


spectral index to be Gaussian priors of a;. The optical depth is sampled by following the tabulated data of median 
value and uncertainties of 7 ê, as described in Ref. Dominguez et al. (2011). 

Our results of Bayesian analysis are shown as follows. At 68% confidence level, we find two sets of posteriors of 
(Aj, ai), i-e., Aj = (2.51 0.46) x 10-4 TeV~*m-?s71, a; = —2.13 +0.11 and A; = (2.14 0.43) x 10-4 TeV ~m ?sT}, 
a; = —1.87+ 0.04. We further depict the one-dimensional posterior PDFs of A; in Fig. 1. In the following, we do 


simulations via sampling A; and a; following their posterior PDFs and 7 following the aforementioned tabulated data. 


3. PROBABILITY OF DETECTING TEV PHOTONS 


Based on the above results, we will estimate the probability of predicting that LHAASO observes at least one photons 
~18 TeV from GRB 221009A. During an observation of T = 2000 seconds, the event number of photons in the energy 
range from (1 — Ag/2)E to (1+ Ag/2)E is given by 


N(E) = TF, (E)Ser(E, 0)Ag (E)E , (5) 


where Ag(E) is the energy resolution of LHAASO-WCDA, as provided in Fig. 26 of Chapter 1 in Ref. Cao et al. 
(2019). We consider the energy range from 0.5 TeV to 30 TeV, with an emphasis on 18 TeV. For each given energy 
E, we perform a set of 10° Monte-Carlo simulations. We count the number of models that predict N(E) > 1 and 
compute the corresponding probability via dividing this number by 10°. 

Our results of Monte Carlo simulations are shown in Fig. 2. We find that the probabilities of predicting that 
LHAASO is capable to detect at least one photons of 18 TeV from GRB 221009A, are about 5.2% and 4.1% in the 
case of priors a; = —1.87 + 0.04 and —2.12+ 0.11, respectively, respectively. Amongst the above two spectral indices, 


the former was obtained during 200-800 seconds after trigger, while the latter was obtained during 500-3500 seconds. 
For comparison, the LHAASO data was accumulated within 2000 seconds after trigger. Based on Ref. Nava (2021), we 
notice that the flux of photons in the TeV band significantly decreases with time. If the events observed by LHAASO 
also preserve this property, though such a guess should be tested with the LHAASO data if possible, we expect that 
the majority of events arrived within the first 1000 seconds. We find that in either case the standard physics is still 
capable to interpret ~18 TeV photons from GRB 221009A. 


4. SUMMARY 


In this work, we have investigated the survival rate of very-high-energy gamma rays within the energy range of 
LHAASO, by taking into account the effect of EBL attenuation. In the framework of standard physics, we simulated 


6 http://side.iaa.es/EBL/ 


94 


95 


202210.00103v1 


chinaXiv 


4 


the probability of detecting the ~18 TeV signal from GRB 221009A to be a few percents. Hints of new physics may be 
displayed at 2.580 level if photons above 26 TeV from GRB 221009A are observed by LHAASO. The above conclusions 
might be altered if we consider other measurements of EBL and optical depth, that are beyond the scope of this paper. 
We would leave such a detailed study to future works. Our research method can be straightforwardly generalized to 
study more data sets of LHAASO and other experiments in the future. 
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